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, We discuss the empirical evidence for a universal Color Glass Condensate and outline 

prospects for further studies at future colliders. Some ramifications for initial conditions 
in heavy ion collisions are pointed out. 

^ ! 1 Introduction 

QCD has been called the perfect theory [T]; as a renormalizable field theory whose validity 
• could extend up to the grand unification scale, it provides the mechanism for generating 

nearly all the mass of the visible universe. The current quark masses are the only external 
Q^' parameters in the theory. Quenched QCD, without dynamical quarks, explains the hadron 

D I spectrum to an accuracy of 10%. These lattice results suggest that gluons play a central 

role in the structure of matter. 

The role of glue in QCD is best understood in the asymptotic weak coupling regimes of 
CsJ ■ the theory where analytical computations are feasible. Much of the discussion in perturbative 

\ QCD (pQCD) has been in the Bjorken-Feynman asymptotics where — *■ oo, s — > oo 

■ and = Q'^/s = fixed. The machinery of precision physics in QCD such as the operator 
QQ I product expansion and factorization theorems are derived in this limit of the theory. The 

progress in this direction has been truly remarkable [3]. In DIS for instance, both coefficient 
functions and splitting functions have been derived to next-to-next-to leading order (NNLO). 

■ What does the hadron look like in the Bjorken-Feynman asymptotics? The DGLAP 
I evolution equations tell us that the gluon distribution grows rapidly with increasing at 
' small XBy However, the phase space density (in a particular gauge and frame), decreases 

rapidly with increasing Q^. The proton become more "dilute" even though the number of 
partons increases; the typical size of resolved partons decreases as faster than the 

increase in the number through QCD evolution. The more dilute the hadron, the cleaner 

■ will be the QCD background for new physics beyond the standard model. 
Much of the current focus in QCD is in quantifying this background. It would be un- 
fortunate however if this were the only focus in QCD studies because the theory, even in 
the weak coupling domain, contains rich and non-trivial dynamics. We speak here of the 
Regge-Gribov asymptotics where xbj — > 0, s — > oo and = fixed. This regime of strong 
color fields is responsible for the bulk of multiparticle production in QCD. What does the 
hadron look like in the Regge-Gribov asymptotics ? The BFKL equation, which resums the 
leading logarithms in x, indicates that the gluon distributions grow even more rapidly in 
this asymptotics. Unlike the Bjorken-Feynman case, the phase space density in the hadron 
grows rapidly as well. The stability of the theory requires that the phase space densities 
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(or more generally, the field strengths squared) be no larger than ~ 1/as- In the pQCD 
framework, mechanisms for the saturation of the growth in the phase space density are 
provided by "higher twist" recombination and screening contributions |4j. These counter 
the bremsstrahlung growth of soft gluons described by the DGLAP and BFKL equations. 
The saturation scale Qs{x) generated by the dynamics demarcates the separation between 
the linear and non- linear regimes of the theory: for momenta Q^, non-linear QCD 

dynamics is dominant, for momenta ^ Ql, weak coupling physics is governed by the 
DGLAP/BFKL evolution equations. 

The universal properties of gluons in the non-linear regime are described by a classical 
effective field theory of dynamical gluon fields coupled to static, stochastic sources. This 
is the Color Glass Condensate (CGC) [S]. The evolution of multi-parton correlators with 
energy is described by the Wilsonian JIMWLK renormalization group (RG) equations [6] . In 
the limit of large nuclei and large Nc, one recovers the Balitsky-Kovchegov equation for the 
forward dipole cross-section [7]. A universal saturation scale arises naturally in the theory 
and its energy dependence is given by the JIMWLK/BK equations. The typical momentum 
of gluons ^ Qs ^ Aqcd ; the bulk of the contributions to high energy cross-sections may be 
therefore described in a weak coupling framework. 

The saturation scale also grows with the nuclear size. A fast compact probe of size 
1/Qs < Rp, where Rp is the proton size, will interact coherently at high energies with 
partons localized in nucleons all along the nuclear diameter. The field strength squared 
experienced by the probe is therefore enhanced parametrically by a factor proportional to 
the nuclear diameter ^ A^^^. As clearly illustrated in the CGC effective theory, the dynamics 
of partons at small x is universal regardless of one speaks of hadrons or nuclei; the latter, 
as we will discuss further, are therefore an efficient (and cheaper) amplifier of the non-linear 
dynamics of these gluons. 

This talk is organized as follows. We will outline our current (limited) understanding of 
the different dynamical regimes in high energy QCD from experiments at HERA and RHIC. 
We will then discuss how experiments at the LHC and future DIS experiments on nuclei can 
help further quantify our understanding. Finally, to illustrate the scope of these studies, 
we will discuss how the strong color field dynamics of partons in nuclear wavefunctions 
contributes to a quantitative understanding of the formation and subsequent thermalization 
of a strongly interacting "glasma" in heavy ion collisions. 

2 The evidence for the CGC from e+p DIS 

A strong hint that semi-hard scales may play a role in small x dynamics at HERA came 
from "geometrical scaling" of the HERA data 0. The inclusive virtual photon-f proton 
cross-section for x < 0.01 and all available scales as a function of t = Q^/Ql, where 
Ql{x) = exp(Ar) GeV^. Here Y = ln(a;o/a;) is the rapidity; xq = 3 ■ lO""* and A = 0.288 
are parameters fit to the data [3 [9]. Further, the inclusive diffractive, vector meson and 
DVCS cross-sections at HERA, with a slight modification [f|in the definition of r, also appear 
to show geometrical scaling [9]. Geometrical scaling of the c4-p data is shown in Fig. [TJ 
A recent "quality factor" statistical analysis [10] indicates that this scaling is robust; it is 
however unable to distinguish between the above fixed coupling energy dependence of Qg 

''The E665 data are a notable exception. 

'^Tjj^VM = iQ'^ + ^'^)/Qbj where M denotes the mass of the diffractive/ vector meson final state. 
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and the running coupling Qs{x) oc cxp(\/y) dependence of the saturation scale. Geometrical 




Figure 1: Geometrical scaling fully inclusive, diffr active and exclusive vector meson cross- 
sections. From 

scaling is only asymptotic in both fixed and running coupling evolution equations 0. Pre- 
asymptotic corrections have been computed previously, to good approximation, in both 
fixed and running coupling cases for the BK equation [TT]. A recent NLO BK analysis [T2] 
suggests that the onset of the scaling asymptotics may be precocious 0, thereby providing a 
possible explanation for its manifestation in the HERA data. A caveat that has been raised 
is that there is a strong correlation between x and in the HERA data [T3] . The scaling 
however persists even where there is a significant lever arm in for small x. Nevertheless, 
geometrical scaling alone is not sufficient evidence of saturation effects and it is important 
to look at the data in greater detail in saturation/CGC models. 

All saturation models [15| express the inclusive virtual photon+proton cross-section as 
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Here 



represents the probability for a virtual photon to produce a quark- 
anti-quark pair of size r — \r±\ and -^2^i^±, x,h±) denotes the dipole cross section for 
this pair to scatter off the target at an impact parameter b^. The former is well known 
from QED, while the latter represents the dynamics of QCD scattering at small x. A simple 
saturation model (known as the GBW model [T^ ) of the dipole cross section, parametrized 

as = 2(l-e"'''<3B,p(^)/4) where Qlpix) = (xq/x)^ GeV^ gives a good qualitative fit to 

the HERA inclusive cross section data for xq — 3 ■ 10~^ and A ~ 0.288. Though this model 
captures the qualitative features of saturation, it does not contain the bremsstrahlung limit 
of perturbative QCD (pQCD) that apphes to small dipoles of size r ^ 1/Qs{x). 



"^The effect of "pomeron loops" on this scaling will be discussed later. 
''For another very interesting take on precocious scaling, see Ref. |14| . 
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In the classical effective theory of the CGC, one can derive, to leading logarithmic ac- 
curacy, the dipole cross section |17| containing the right small r limit. This dipole cross 
section can be represented as [T8] 

d ^ 

^ = 2 [1 - exp (-r2F(x, r)T,{h^))] , (2) 

where Tp{h±) is the impact parameter profile function in the proton, normalized as / d^bj^ Tpih 
1 and F is proportional to the gluon distribution [19| 

F{x, r2) = TT^a, + 4/r2) xg {x, ^ll + A/r'') l{2N,) , (3) 

evolved from the initial scale /^o by the DGLAP equations. The dipole cross section in 
Eq. ^ was implemented in the impact parameter saturation model (IPsat) |18| where the 
parameters are fit to reproduce the HERA data on the inclusive structure function i^2- Here 

is defined as the solution of ^^{x,r^ = l/Ql{x,h^_)) = 2(1 - e-i/^)0. 

The IPsat dipole cross section in Eq. ([2]) is valid when leading logarithms in x in pQCD 
are not dominant over leading logs in . At very small x, where logs in x dominate, 
quantum evolution in the CGC describes both the BFKL limit of linear small x evolution as 
well as nonlinear JIMWLK/BK evolution at high parton densities [6l[7j. These asymptotics 
are combined with a more realistic 6-dependence in the b-CGC model PHI EI] • Both the 
IPsat model and the b-CGC model provide excellent fits to HERA data for x < 0.01 [HI [22]. 
An important caveat [23] to the success of the saturation models is that the saturation scale. 




Figure 2: The saturation scale vs l/x in the IPsat and b-CGC models [21]. 

at median impact parameters, extracted from these fits is < 1 GcV^ even at the lowest x 
values at HERA [21] [24] . The saturation scale extracted from the fit in the IPsat model 
is shown in Fig. [2| We should note however that the uncertainty in the magnitude of the 
saturation scale is significant and is a factor of 2 larger in recent CGC fits [55]. NLO 

' This choice of is equivalent to the saturation scale in the GBW model for the case of a Gaussian dipole 
cross section. 
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computations in the small x dipole framework are now becoming available [26] : these will 
provide theoretical guidance into precisely how the coupling runs as a function of Qs at 
small X. Finally, from Fig. [2l it is clear that the energy dependence of the extracted Qs is 
significantly stronger than those predicted in non-perturbative models [27] . 

3 The evidence for the CGC from e+A DIS and d+A and A+A 
collisions 

The strong field dynamics of small x partons is universal and should be manifest in large nu- 
clei at lower energies than in the proton. In Fig. [3] (left), we show the well known shadowing 
of in the fixed target e+A E665 and NMC experiments. Expressed in terms of t = /Ql 
(Fig. [3] (right)), the data show geometrical scaling [28]. In Ref. [28], the A dependence of 
Qs is determined to be A^/^ and not A^^^ as suggested in a simple random walk picture. 
However, as we shall discuss shortly, this conclusion is a little misleading. A study of nuclear 
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Figure 3: Left: Shadowing of F2 from the NMC and E665 fixed target experiments. Right: 
The data scaled as a function of r = Q'^/Q^ [28] . 



DIS in the IPsat CGC framework was performed in Ref. 
dipole cross section is well approximated by ' 
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The average differential 
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, where 



TaO-^±) is the well known Woods Saxon distribution. Here cr^jp is determined from the IPsat 
fits to the e+p data; no additional parameters are introduced for eA collisions. In Fig. [J] 
(left), the model is compared to NMC data on Carbon and Calcium nuclei-the agreement 
is quite good. In Fig. |4] (right), we show the extracted saturation scale in nuclei for both 
central and median impact parameters. To a good approximation the saturation scale in 
nuclei scales as Qsj^{x,brtwd ) ~ Qs.pi^T^mcd.) ■ (A/x)^/"^. The factor of 200^/'^ ~ 6 gives a 



^This is considerably larger than the simplest estimate of a 
in the GBW model, which yields « ~ A^/^- " 
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huge "oomph" in the parton density of a nucleus relative to that of a proton at the same 
X. Indeed, one would require a center of mass energy ~ 14 times larger^ in an e+p collider 
relative to an e+Au collider to obtain the same j^{bmed.)- The reasons for the additional 
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Figure 4: Left: Comparison of the IPsat model (with no adjustable parameters) to the NMC 
data. Right: The A and x dependence of the saturation scale in the IPsat model pS] , 

enhancement are two fold. Firstly, because the density profile in a nucleus is more uniform 
than that of the proton, Qg p(&med.) is only ~ 35% of the value at & = 0; in contrast, in 
gold nuclei it is 70%. Because the median impact parameter dominates inclusive scattering, 
this effect gives a significant enhancement to the effective Qs- The second reason for the 
enhancement is the DGLAP-like growth of the gluon distribution in the IPsat nuclear dipole 
cross section. For two nuclei, A and B (with A > B), in a "smooth nucleus" approxima- 
tion ij:t=iTpihi_ - b^,) ATAih^)), gl^ « F^Sitr "^^""^'^ ^ ^""^ '^"^^'"^ 
eq. [3l The scaling violations in F imply that, as observed in Refs. [18l[3T], the growth of 
Qs is faster than A^^^. Also, because the increase of F with is faster for smaller x, the 
A-dependence of Qs is stronger for higher energies. In contrast, the dipole cross section in 
the b-CGC model depends only on the combination^ rQs{x) without DGLAP scaling viola- 
tions. It therefore does not have this particular nuclear enhancement. Another interesting 
possibility, following from running coupling corrections to the leading logs in x, is that QCD 
evolution actually depletes the nuclear enhancement of Qs at very small x [30|. Precise 
extraction of the A dependence of Qs can therefore help distinguish between "classical" and 
"quantum" RG evolution at small x. 

We now turn to a discussion of CGC effects in hadronic collisions. A systematic treatment 
of the scattering of two strong color sources (such as two high energy nuclei) is discussed in 
Sectional To leading order, the problem reduces to the solution of the classical Yang-Mills 
(GYM) equations averaged over color sources for each nucleus [211 [33]; the variance of this 
distribution of sources is proportional to Qg a- Besides the nuclear radius, Q^ ^ the only 
scale in the problem, and the Qg a '5s,p ' {A/x)^-^ expression for the saturation scale 
was used in CGC models of nuclear collisions to successfully predict the multiplicity [33] 
and the centrality dependence of the multiplicity [34] dependence in gold-|-gold collisions at 

''At extremely high energies, this statement must bo qualified to account for the cfl'ccts of QCD evolu- 
tion [30l . 

'With the caveat that it has BFKL-like violations that vanish asymptotically with Y . 
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RHIC. The universality of the saturation scale also has a bearing on the hydrodynamics of 
the Quark Gluon Plasma (QGP); the universal form leads to a lower eccentricity [37| (and 
therefore lower viscosity) than a non-universal form that generates a larger eccentricity j38j 
(leaving room for a larger viscosity) of the QGP. 

For asymmetric (off-central rapidity) nuclear collisions, or proton/deuteron-|- heavy nu- 
cleus collisions, A: j_ -factorization can be derived systematically for gluon production, at lead- 
ing order, in the CGC framework 09]. The simplicity of k± factorization is convenient for 
phenomenology; predictions based on this formalism describe the rapidity distributions in 
A + A collisions [35] and the phenomenon of "limiting fragmentation" ^0]. The latter, 
and deviations thereoff, are described by solutions of the BK-equation. Predictions for the 
multiplicity distribution in A+A collisions at the LHC [41] for both GBW and classical 
CGC (MV) initial conditions | give a charged particle multiplicity of 1000-1400 in central 
lead+lead collisions at the LHC. The results are shown in Fig.[5j In deuteron+gold collisions 




Figure 5: Prediction for limiting fragmentation and deviations away from it at LHC energies. 
The bands denote the range in the predictions for GBW and MV models. From [4T] . 

at RHIC, the normalized ratio RpA of the inclusive hadron spectrum relative to the same 
in proton-|-proton collisions shows a mild "Cronin" peak at mid-rapidities corresponding 
to multiple scattering in the classical CGC; at forward rapidities, however, RpA decreases 
systematically below unity. In the CGC, this reflects quantum evolution of the dipole cross- 
section in a large nucleus and has the same origin as the extension of the geometrical scaling 
regime gg to Q > Qs- This effect 0, should also exist in hadronic collisions [13] ; specifically, 
it was predicted this would occur in deuteron-|-gold collisions [33]. In general, RpA while 
suggestive, is not an ideal variable because it is not clear the same formalism applies to p-l-p 
collisions at the same rapidity. Data on the inclusive hadron spectrum in deuteron-|-gold col- 
lisions can be directly compared to model predictions [45] Q. The result is shown in Fig. (6] 
For a comprehensive review of applications of CGC picture to RHIC phenomenology, we 

JThe McLerran-Vcnugopalan (MV) initial condition has the same form as the IPsat dipole cross-section 
discussed earlier. 

'^Quantum evolution here corresponds to the BK anomalous dimension of 7 = 0.63 in the dipole cross- 
section, as opposed to 7 = 1 (DGLAP) and 7 = 0.5 (BFKL). 

'The same analysis also gives good agreement for the forward p+p spectrum at RHIC |46) . 
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refer the reader to Ref. [17]. There are a couple of caveats to this picture. Firstly, k± fac- 
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Figure 6: The inclusive fc^ distributions in deuteron+gold collisions compared to theory 
curves for different rapidities. From |45j . 

torization is very fragile. It docs not hold for quark production even at leading order in the 
parton density [15], albeit it may be a good approximation for large masses and transverse 
momenta [13]. For gluon production, it does not hold beyond leading order in the parton 
density [5D1 [33] . Secondly, a combined comprehensive analysis of HERA and RHIC data is 
still lacking though there have been first attempts in this direction [5T] . 

4 The future of small x physics at hadron colliders and DIS 

The LHC is the ultimate small x machine in terms of reach in x for large Q^. A plot from 
Ref. [S2] illustrating this reach is shown in Fig. [7] (left). For a recent review of the small 
X opportunities at the LHC, see Ref. [S3]. The LHC will provide further, more extensive 
tests of the hints for the CGC seen at RHIC. At very high energies, a novel "diffusive 
scaling" regime has been proposed, which incorporates the physics of Pomeron loops [54] . 
Recent developments were reviewed at DIS06 by lancu and at DIS07 by Shoshi [55] ; possible 
signatures at the LHC have been proposed |56| . However, very recent computations including 
running coupling effects suggest that this regime is unlikely to be accessed realistic collider 
energies [57] . 

The universality of parton distributions is often taken for granted but factorization the- 
orems proving this universality have been proven only for a limited number of inclusive 
final states. However, as we have discussed, small x is the domain of rich multi-parton 
correlations. These are more sensitive to more exclusive final states for which universality 
is not proven [58j . Therefore, while the LHC will have unprecedented reach in x, precision 
studies of high energy QCD and clean theoretical interpretations of these motivate future 
DIS projects. Two such projects discussed at this conference are the EIC project in the 
United States [59] and the LHeC project in Europe |60] . 

As we discussed previously, strong color fields may be more easily accessible in DIS off 
nuclei relative to the proton due to the "oomph" factor. In Fig. [7] (right), we show the 
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Figure 7: Left: Kinematic x-Q^ reach of different final states at the LHC compared to other 
experiments with nuclei . From |52| . Right: The saturation scale in the proton, calcium 
and gold in the kinematic acceptance of the EIC. 

saturation scale j^{x) overlaid on the x-Q^ kinematic domain spanned by the EIC. It 
is interesting that there is a significant kinematic domain where Qg ^ > Q^, including in 
particular ^ > I GeV^. In the weak field regime where ^ p^, we are accustomed to 
thinking of as = as{Q^)- In the strong field regime, where Qg a ^ Q^, we likely have instead 
as = asiQsA)- suggested by the figure, the EIC (and clearly the LHeC) will cleanly 
probe the cross-over regime from linear to non-linear dynamics in QCD. A particularly 
striking feature of e-f A DIS will be diffractive scattering [611 [29] ; it is anticipated that ~ 30% 
of the cross-section corresponds to hard diffractive final states. For further discussion of the 
physics of an Electron Ion collider, see Ref. [62] , 

5 From CGC to QGP: how classical fields decay in the exploding 
Glasma 

The word "Glasma" describes the strongly interacting matter in heavy ion collisions from 
the time when particles are produced in the shattering of two CGCs to the time when a 
thermalized QGP is formed [64] ■ We will discuss here a systematic approach to computing 
particle production in heavy ion collisions to NLO. This approach suggests a deep connection 
between quantum evolution effects in the nuclear wavcfunction and instabilities that may 
be responsible for fast thermalization of the Glasma. 

A cartoon of multi-particle production in a heavy ion collision is shown in Fig. [8l The 
probability of producing n particles, in field theories (such as the CGC) with strong external 
sources can be expressed as [55] 

K^e.J-l,Y.^.)±l ■£ '-^^ (4) 

\ ^ r / p=l ^ Cti + ---Qp=Tl ^ 

where br denotes the sum of vacuum-to- vacuum graphs with r cuts. This formula has 
remarkable features: a) P„ is non-perturbative in g even for g <C 1-no simple power expansion 
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Figure 8: Cartoon of gluon production in the collision of two sheets of Colored Glass. The 
dots denote large x color sources. 

in terms of g exists, b) P„, for any n, gets contributions from cut tree vacuum graphs-this 
would not apply for field theories in the vacuum, c) Even at tree level, P„ is not a Poisson 
distribution which is counter to presumptions that classical field theories have only "trivial" 
Poissonian correlations. The simple formula in Eq. [4] contains many features of the well 
known AGK calculus [66] of multi-particle production. 

Computing the probabilities in Eq.|3]is hopeless even for g <^1. Fortunately, a systematic 
expansion in powers of g exists for moments of the multiplicity. Both the LO and NLO 
multiplicity can be represented in terms of solutions of equations of motion with retarded 
boundary conditions. At leading order, these are solutions A^^^ to the Yang-Mills equations; 
these equations, with boost invariant CGC initial conditions, were solved numerically in 
Ref. [321 to compute the gauge fields at late times. 

A next-to-leading order (NLO) computation is important to understand the renormaliza- 
tion and factorization issues that arc fundamental to any quantum field theory. As we shall 
discuss shortly, it is also important to understand the quantum fiuctuations that generate 
the plasma instabilities which may speed up thermalization. Remarkably, the NLO contri- 
butions can be computed by solving the initial value problem of small fiuctuation equations 
of motion with retarded boundary conditions. A similar algorithm has been constructed and 
implemented to study quark pair production in the classical CGC background field I65j. 

In the Glasma, the classical LO boost invariant E and B fields are purely longitudinal 
at T = 0. The corresponding momentum distributions, at r > \/Qs, are very unstable- 
indeed, they lead to an instability which may be analogous to the well known Weibel in- 
stability in electromagnetic plasmas. For a review and relevant references, see Ref. |67j . 
3-l-l-D numerical simulations demonstrate that small rapidity dependent quantum fluctua- 
tions grow exponentially and generate longitudinal pressure |68| . The initial "seed for the 
simulations corresponded to "white noise" Gaussian random fluctuations. The maximally 
unstable modes of the longitudinal pressure grow as cxp (C-y/Avr) with C ~ 0.425; this form 
of the growth was previously predicted for Weibel instabilities in expanding plasmas |69j . 
Albeit the solutions of the Yang-Mills equations display similar features to the HTL studies, 
a deeper understanding of this connection is elusive. 

First quantum corrections to the classical background field of two nuclei at r = give 
initial conditions j7D] that are quite different from those in Ref. |68| . Simulations are under- 
way to determine whether these initial conditions speed up thermalization. A full treatment 
of quantum fiuctuations requires that we understand how some NLO contributions are ab- 
sorbed in the evolution of the nuclear wavefunctions with energy while the rest contribute 
to gluon production. A proof of this high energy "factorization" is in progress |71j . To 
fully understand fast thermalization in the presence of instabilities, one also needs a kinetic 
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theory of the Glasma that describes the decay of classical fields into particles. A first step 
has been made in this direction [72] . 
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